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An extended Keldysh formalism, well suited to properly take into account the initial correlations, 
is used in order to deal with the time-dependent current response of a resonant tunneling system. We 
use a partition-free approach by Cini in which the whole system is in equilibrium before an external 
bias is switched on. No fictitious partitions are used. Despite a more involved formulation, this 
partition-free approach has many appealing features being much closer to what is experimentally 
done. In particular, besides the steady-state responses one can also calculate physical dynamical 
responses. In the noninteracting case we clarify under what circumstances a steady-state current 
develops and compare our result with the one obtained in the partitioned scheme. We prove a 
Theorem of asymptotic Equivalence between the two schemes for arbitrary time-dependent distur- 
bances. We also show that the steady-state current is independent of the history of the external 
perturbation (Memory Loss Theorem). In the so called wide-band limit an analytic result for the 
time-dependent current is obtained. In the interacting case we work out the lesser Green function 
in terms of the self energy and we recover a well known result in the long-time limit. In order to 
overcome the complications arising from a self energy which is nonlocal in time we propose an exact 
non- equilibrium Green function approach based on Time Dependent Density Functional Theory. 
The equations are no more difficult than an ordinary Mean Field treatment. We show how the 
scattering-state scheme by Lang follows from our formulation. An exact formula for the steady- 
state current of an arbitrary interacting resonant tunneling system is obtained. As an example the 
time- dependent current response is calculated in the Random Phase Approximation. 

PACS numbers: 05.60.Gg Quantum transport 
72.10.Bg General formulation of transport theory 

85.30.Mn Junction breakdown and tunneling devices (including resonance tunneling devices) 



I. INTRODUCTION 

A resonant tunneling system is essentially a mesoscopic 
region, typically a semiconductor heterostructure, cou- 
pled to two metallic leads, which play the role of charge 
reservoirs. In a real experiment the whole system is in 
thermodynamic equilibrium before the external distur- 
bance is switched on and one can assign a unique tem- 
perature and chemical potential /i. Therefore, the 
initial density matrix is p ~ exp[—P{H — fj,N)] where H is 
the total Hamiltonian and N is the total number of par- 
ticles. By applying a bias to the leads at a given time, 
charged particles will start to flow through the central 
device from one lead to the other. As far as the leads are 
treated as noninteracting, it is not obvious that in the 
long-time limit a steady-state current can ever develop. 
The reason behind the uncertainty is that the bias rep- 
resents a large perturbation and, in the absence of dissi- 
pative effects, e.g., electron-electron or electron-phonon 
scatterings, the return of time-translational invariance is 
not granted. 

An alternative approach to this quantum transport 
problem has been suggested by Caroli et al.^-^ who state: 
"It is usually considered that a description of the sys- 
tem as a whole does not permit the calculation of the 
current" i Their approach is based on a fictitious par- 
tition where the left and right leads are treated as two 
isolated subsystems in the remote past. Then, one can fix 



a chemical potential /Iq and a temperature (3^ for each 
lead, a — L, R. In this picture the initial density matrix 
is given by p exp[-/3L(i7L - fiLNL)]exp[~l3R{Hii - 
fifiNfj)], where Hlr and N^^n now refer to the iso- 
lated L, R lead. The current will flow through the sys- 
tem once the contacts between the device and the leads 
have been established. Hence, the time-dependent per- 
turbation is a lead-device hopping rather than a local 
one-particle level-shift. Since the device is a mesoscopic 
object, it is reasonable to assume that the hopping per- 
turbation does not alter the thermal equilibrium of the 
left and right charge reservoir and that a non-equilibrium 
steady state will eventually be reached. This argument 
is very strong and remains valid even for noninteracting 
leads. Actually, the partitioned approach by Caroli et 
al. was originally applied to a tight-binding model^ de- 
scribing a metal-insulator-metal tunneling junction and 
then extended to the case of free electrons subjected 
to an arbitrary one-body potential.^ This extension was 
questioned by Feuchtwang)^ the controversy was about 
the appropriate choice of boundary conditions for the 
uncontacted-system Green functions. In later years the 
non-equilibrium Green function techniques''^ in the par- 
titioned approach framework were mainly applied to in- 
vestigate steady-state situations. An important break- 
through in time-dependent non-equilibrium transport 
was achieved by Wingreen et aZ.^'^'^'^° Still in the frame- 
work of the partitioned approach, they derive an expres- 
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sion for the fully nonlinear, time-dependent current in 
terms of the Green functions of the mesoscopic region 
(this embedding procedure holds only for noninteracting 
leads). Under the physical assumption that the initial 
correlations are washed out in the long-time limit, their 
formula is well suited to study the response to external 
time-dependent voltages and contacts. 

The limitations of the partitioned approach are essen- 
tially three. First, it is difficult to partition the electron- 
electron interactions between the leads and between the 
leads and the device. These interactions are important 
for establishing dipole layers and charge transfers which 
shape the potential landscape in the device region. Sec- 
ond, there is a crucial assumption of equivalence between 
the long-time behavior of the 1) initially partitioned and 
biased system once the coupling between the subsystems 
is established and 2) the whole partition-free system in 
thermal equilibrium once the bias is established. Third, 
the transient current has no direct physical interpretation 
since in a real experiment one switches on the bias and 
not the contacts; moreover, there is no well defined pre- 
scription which fixes the initial equilibrium distribution 
of the isolated central device. 

In this paper we use a partition-free scheme without 
the above limitations. This conceptually different time- 
dependent approach has been proposed by Ciniiii He 
developed the general theory for the case of free elec- 
trons described in terms of a discrete set of states and 
a continuum set of states with focus on semiconductor 
junction devices. For a one-dimensional free-electron 
system subjected to a time-dependent perturbation of 
the form UQ{t)<d{x), where U is the applied bias and 
{x, t) is the space-time variable, the Cini theory yields a 
current-voltage characteristics which agrees with the one 
obtained by Feuchtwang2»^ in the partitioned approach. 
This result is particularly important since it shows that 
a steady state in a partition-free scheme develops even 
in the noninteracting case. Moreover, it demonstrates 
an equivalence which had previously been assumed. In 
the present work we extend the partition-free approach 
to noninteracting resonant tunneling systems and also 
to interacting such systems - in both cases using arbi- 
trary time-dependent disturbances. We shall clarify un- 
der what circumstances a non-equilibrium steady state 
can develop and discuss the equivalence of the current- 
voltage characteristics obtained by Jauho et al~ and that 
obtained by us. One of the advantages of the partition- 
free scheme over the traditional methods lies in the ability 
of the former to calculate transient physical (i.e. mea- 
surable) current responses. 

The plan of the paper is the following. In Section 
im we develop the general formalism which properly ac- 
counts for the initial correlations. We derive a solution 
of the Keldysh equations for the lesser and the greater 
Green functions in noninteracting and interacting sys- 
tems. An exact and alternative treatment based on 
Time Dependent Density Functional Theory^^ (TDDFT) 
is proposed in order to calculate the total nonlinear time- 



dependent current. The current response of a noninter- 
acting resonant tunneling system is discussed in Section 
mil We specify when the partitioned and the partition- 
free schemes yield the same asymptotic current (Theo- 
rem of Equivalence) and how this current may depend 
on history (Memory Loss Theorem). The general re- 
sults are illustrated by model calculations. In Section 
II VI we consider an interacting resonant tunneling system 
with interacting leads. The TDDFT approach is com- 
pared with earlier works by Lang et alM^i^ and Taylor et 
ali^^ Assuming that a steady state is reached we write 
down an exact formula for the nonlinear steady-state cur- 
rent. As a simple example we also study the current re- 
sponse in the Random Phase Approximation (RPA) of 
a capacitor-device-capacitor junction. Our main conclu- 
sions are summarized in Section Ivl 

II. GENERAL FORMULATION 

A. Noninteracting Systems in the Presence of an 
External Disturbance 

Let us consider a system of noninteracting electrons 
described by an unperturbed Hamiltonian 

^ ^ -^m.n^^^n? i_'-^^m,n — T^m.n (1) 

mn 

and by a time-dependent disturbance of the form 

Hu{t) ^^UmA^)cLcn, [U{t)]m,n^UmA^)^ (2) 
ran 

with U{t) = for any t < to. In Eqs. C3l-©i c„i, cjj are 
Fermi operators in some suitable basis, and we use bold- 
face to indicate matrices in one-electron labels. Without 
loss of generality one can take to = 0. The system is in 
equilibrium for negative times. 

1. Elementary Derivation 

We first obtain the Green function by elementary 
means without resorting to any Keldysh techniques. For 
a noninteracting system everything is known once we 
know how to propagate the one-electron orbitals in time 
and how they are populated before the system is per- 
turbed. The time evolution is fully described by the re- 
tarded or advanced Green functions G^'^, and the initial 
population at zero time, i.e., by G^(0; 0). The real-time 
Green functions are defined by 

G^-^it; t') = ±e(±t T t') [G> [t- 1') - G<{t- 1')] , 

with 

G>At;t') - ~z{cAt)clit')), 
G<At;t') = ^{cUt')cAt)), 
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where the operators are Heisenberg operators and where 
the averages are with respect to the equihbrium grand- 
canonical ensemble. Because there are no inter-particle 
interactions, the equation of motion for the electron op- 
erators simplifies to 

n 

where K{t) = T + U{t) is the full one-body Hamilto- 
nian matrix. Consequently, the time evolution of c„j is 
given by the one-electron evolution matrix S{t), Cm{t) = 
E„ 5'm,n(i)c„(0), where S obeys iS{t) = K{t)S{t), with 
initial value <S'(0) = 1. We insert the time-evolved oper- 
ators in the definitions of the G matrices to obtain 

G^'^it; t') = T«e(±t T t')S{t)SHt'), (3) 

and 

G^(t;t') = S{t)G'^{0;0)S^t') 

= G^{t;0)G^{0;0)G^{0;t'), (4) 

where the last equality holds for any t,t' > 0. We observe 
that the instantaneous current can be expressed in terms 
of G"^ (t; t) , and thus the problem of finding the current 
is reduced to that of finding the retarded Green function 
and the equilibrium population of the one-electron levels. 
We note in passing that the initial populations can be 
expressed as G^(0;0) = if{T), where / is the Fermi 
fimction. Because T is a matrix, so is f{T). 

The above solution for the lesser /greater Green func- 
tion was also derived by CiniAi with an equation-of- 
motion approach. He also pointed out that they can be 
derived in the framework of the Kcldysh formalismS. as a 
finite-temperature extension of a treatment by Blandin 
et aliL 



2. Derivation Based on the Keldysh Technique 

In this Subsection we give an alternative derivation of 
Eq. I0J using an extension of the Keldysh formalism. 
There are two reasons for giving another derivation. On 
one hand, we will use the Keldysh formalism taking due 
account of the prescribed integration along the imaginary 
axis. This will allow us to understand what kind of ap- 
proximations are made in the partitioned approach. On 
the other hand, the derivation below clearly shows how 
the electron-electron interaction can be included. 

We introduce the Green function 

Gyn,n{z\z') = -i{T[Cm{z)Cn{z')]) (5) 

which is path-ordered on the oriented contour 7 of Fig. 
^ In Eq. (jsj z = t + T is the complex variable running on 
7 with t = 3l[z], T — i'^[z\ while A and B are the end- 
points of 7. Further, Cm{z) and c„(z) are Heisenberg 
operators defined by the non-unitary evolution operator 



for complex times z. They are in general not Hermitian 
conjugates of one another, but the usual equal-time an- 
ticommutation relations {0^(2:), c„(z)} — 5„i^n are still 
obeyed. As before the average is the grand-canonical av- 
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FIG. 1: Contour suited to include the effect of the initial correla- 
tions, see also Section Iff Bl 



erage. On the vertical track going from to — i/3 we have 
K{t) = K{0) = T independent of r. Therefore, the 
Green function satisfies the relations 

G(A;z') = -e'3MG(i?;z') 
G{z-A)^-e-^f'G{z-B)- ^> 

Next, we write the total Hamiltonian H{t) = Hq + 
Hu (t) as the sum of a diagonal term and an off-diagonal 
one 

m mn 

The quantities em{T) = £m and Kn,n(T) = Vm,n are 
constants on the vertical track. [The decomposition 
above is completely general. In our model examples dis- 
cussed later, the diagonal part [£{z)]m,n = Sm,n£miz) 
will represent an uncontacted system and the off-diagonal 
[^(2;)]m,n = Vm,n{z) the contacts.] The Green function 
G{z; z') is obtained by solving the equation of motion 

[and its adjoint] with boundary conditions We define 
g(z; z') as the uncontacted Green function. The g satis- 
fies Eq. (T)) with V = and obeys the same boundary 
conditions of the contacted G. The unique g resulting 
from such a scheme belongs to the Keldysh spaceii and 
can be written as 

g(z; z') = e(z, z')g>{z; z') + e(z', z)g<(z; z'), 

where 0(z, z') = 1 if z is later than z' on 7 and zero 
otherwise. g^(z;z') is analytic for any z later than z' 
while g^(z;z') is analytic for any z' later than z; they 
are given by 

g<(z;z') - if{£)e-^S^"'-^^^-^\ 

e{z;z') = z[/(5)-l]e-'/.''^^-^(^-), (8) 
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where £ = £(0) and the integral appearing in the ex- 
ponential function is a contour integral along 7 going 
from z' to z. Choosing z and z' on the real axis and 
reduce to the real-time lesser and greater component. 
From Eqs. ((SJ one can easily verify that the correspond- 
ing retarded and advanced component can be written as 



g«'^(i;0 = T*e(±iTOe- 



■11 die it) 



(9) 



The uncontacted g allows to convert Eqs. lO into an 
integral equation which preserves the relations ®: 

G(z; z') = g(z; z') + f dz g(z; z)V{z)Giz; z'). (10) 

Using the Langreth theorem.^^ one finds 

5 + ■ v\ -g^ + G^-F-g^ + G^y^gr, (11) 



G^ = 



where we have used the short hand notation • to denote 
integrals along the real axis, going from to cxd, and * for 
integrals along the imaginary vertical track, going from 
to — For the sake of clarity we have also introduced 
the symbols ] and [: any function with the superscript ] 
is intended to have a real first argument and an imaginary 
second argument; the opposite is specified by [. In Eq. 
(inil, z') = S{z - z')V{z); for V we don't need to 
say more since it is always foregone and followed by • or 
* so that no ambiguity arises. In particular we note that 
■kV-k implies a simple matrix multiplication since along 
the vertical track V is a constant matrix times the delta 
function. 

The equation for G> contains G{t; r) with one real and 
one imaginary argument. This coupling does not allow to 
get a closed equation for G{t; t') with two real arguments, 
unless V = on the vertical track. Conversely, G^ and 
G^ satisfy an integral equation without any coupling: 



G = g 



R.A 



G 



R,A 



y .g 



R,A 



(12) 



Eq. (|11|) can be solved for G^ and one obtains 



G^ 



5 + G'^ ■ V 



6 + V -G' 



-G'^ i^Vi^s} ■\6 + V -G 



(13) 



From Eq. JSJ and Eq. © we have g^{t;t') = 
gi^(t;0)g^(0;0)gA(0;t') and g(T;t) = -zg(T; 0)gA(0; t), 
so that Eq. (|13|l can be rewritten as 

G^{t;t') = G«(t;0)g^(0;0)G^(0;i') 

g1 ★y^gl (<;0)G^(0;t')- (14) 



The above expression for G> deserves a brief comment. 
Indeed, the first term on the r.h.s. is exactly what 
one got in the partitioned approach, where the hopping 
parameters Vm.n vanish along the vertical track. It is 



usually argued that if t,t' 00 the second term van- 
ishes. However, we point out that in the noninteract- 
ing case this is not true. If in the long-time limit some 
physical response functions, e.g., the current, are cor- 
rectly reproduced by using the partitioned G^{t;t') — 
G"(t;0)g^(0;0)G'^(0;i') other kind of argumentations 
should be invoked. We shall come to this point later on. 

To proceed further we need the Dyson equation for 
G{t; t). Exploiting the identity g(i; r) = zg^(t; 0)g(0; r), 
we find 

G{t- t) = iG^'it- 0)g(0; r) + [gU y ★ g] (i; r). (15) 

Eq. (|15|l can be solved for G{t\ r). From the Dyson equa- 
tion G(t;t') = g(r;T') + [G * F * g] (r; r'), it follows 
that [5 - y * g]"^ (f; r) = [(5 + F * G] (f; r) and hence 

G{t-T)^iG^{t-Q)G{Q]T). (16) 

Substituting Eq. (fTE|l into Eq. ifTHl one gets 

G^{t- 1') = G^{t- 0)G^(0; 0)G^(0; t'). (17) 

Eq. ((T7|l coincides with Eq. as it should. 



B. Interacting Systems in the Presence of an 
External Disturbance 

In the interacting case we keep track of the interac- 
tions by introducing a self-energy matrix. Then, Eq. ^ 
becomes 



i--£{z)-V{z)-i:\z)^G{z;z') 
^5[z-z')+ I dzT,c{z;z)G{z-z'). 



(18) 



Here S is the self-energy part which is local in time 
and it consists of a Hartree and an exchange term. The 
remaining part of the self energy contains the contri- 
butions coming from the correlation and belongs to the 
Keldysh spaceti^ 

Sc(z; z') = 6(2, 2')S> (2; z') + e(z', 2)S<(z; z'). 

Like G, the self energy and its components are matrices 
in the one-electron labels. No simple expressions, like 
©"O, can now be directly obtained from the equation 
of motion and the Keldysh formalism is unavoidable. 

A proper treatment of the initial correlations naturally 
leads to an extension of the Keldysh equations. The gen- 
eralization was put forth by Wagner=22i who obtained a 
minimal set of five independent integro-differential equa- 
tions for the unknowns G^, G^, G< (or G>) (or the 
Keldysh Green function G^ = G> + G<), G^ (or G^) 
and the thermal Green function G with two imaginary 
arguments. In Appendix fXl we exploit the results of the 
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previous Subsection to prove that the solution for 
can be written as 

G^(<; t') = G^{t- 0)G^(0; 0)G^(0; t') + A^{t; t'), (19) 

where 

A^{t;t') = tG^{t-0)G>iO-t')^tG<{t;0)G^{0;t') 
- G^(i;0)G^(0;0)G^(0;i') 



G^ • [s^ + Sl *G*Sr 



G' 



(t-X). 



This resuh clearly reduces to Eq. (|17|l if the self energy 
vanishes since G^(0;0) = [G^(0;0)]t = -i. We observe 
that if the Green functions vanish when the separation 
of their time arguments goes to infinity, Eq. I|19|l yields 
a well known identity 



lim G^{t;t') 



G' 



S^G^ 



(20) 



Eq. I|2U|I is well suited to study the long-time response 
of an interacting system subjected to an external time- 
dependent disturbance. On the other hand, if one is in- 
terested in the short-time response Eq. ifTUIl cannot be 
simplified. In some cases it might be simpler to use an 
alternative approach. Below we propose an exact non- 
equilibrium Green function treatment based on TDDFT 
and discuss the relations to ordinary Mean Field approx- 
imations. 



C. Mean Field Theory and Relations to TDDFT 

Any Mean Field Theory is a one-particle-like approx- 
imation in which each particle moves in an effective av- 
erage potential independently of all other particles. The 
mean-field potential is local in time, meaning that Sc is 
discarded. Consequently, all the results of the Section 
III Al can be reused provided we substitute K hy K . 
Thus, no extra complications arise if we treat an interact- 
ing system at the Hartree-Fock level. To be specific, let 
us focus on the Coulomb interaction and on paramagnetic 
systems (so that the self energy and the Green function 
are diagonal in the spin indices). Then, it is natural to 
choose the one-electron index as the coordinate r of the 
particle and to split the self energy ^, (z) = S''(r, r', z) 
as a sum of the Hartree and the exchange term 



S'(r,r' 



z) = yH(r,z)(5(r-r') + Sx(r,r',z). 



For extended systems, the Hartree potential Vh and the 
Coulomb potential from the nuclei are separately in- 
finite but with a finite sum. Together with the exter- 
nal field U these terms form the classical electrostatic 
potential Uc — U -\- Va + Vn- The Green function 
Gr.r' {z; z') = G(r, z; r', z') can be obtained from the self- 
consistent solution of the equation of motion and the 



lesser/greater component can be written as 
G|F(r,i;r',t') = J dfdf' G^p{r,t;r,0) 

xG|F(f,0;f',0)G^F(f',0;r',i'), 

where the subscript MF has been used to stress that it is 
a Mean Field approximate result. In the ordinary Many 
Body Theory one has to abandon the one-particle pic- 
ture in order to improve the approximation beyond the 
Hartree-Fock level. This leads to a self energy nonlocal 
in time and hence to the complicated solution (|19|l . 

In the case we only ask for the density 
n{r,t) — ~iG^{r,t;r,t) the original Density 
Functional TheorySLSS and its finite-temperature 
generalization^^ has been extended to time-dependent 
phenomenaiiSiSi The theory applies only to those cases 
where the external disturbance is local in space, i.e., 
Ur.T'{t) — (5(r — r')C/(r, i). For t > we switch on an 
external potential U{r,t) to obtain a density n{r,t). 
The Runge-Gross theorem states that if we instead 
had switched on a different U'{r,t) [giving a different 
n'(r, i)], then n{r,t) = n'{r,t) implies U{r,t) = [/'(r,t). 
Thus U{r,t) is a unique functional of n{r,t). Runge 
and Gross also show that one can compute 'n{r,t) in a 
one-particle manner using an effective potential 

C/°f^(r,<) = ;7c(r,t)+«,c(r,i). 

Here, Uxc accounts for exchange and correlations and is 
obtained from an exchange-correlation action functional, 
Wxc(r, i) = SA^c[n]/Sn{r,t). In our earher language this 
corresponds to an effective self energy which is local in 
both space and time. The TDDFT one-particle scheme 
corresponds to a fictitious Green function Q{r,z;r',z') 
which satisfies the equations of motion Q with [fr,r' + 
Vryiz)] replaced by ^(r - r')hV?/2 + U^^ir,/)]. As a 
consequence we have 

g^{r,t;r',t') = J dfdf' g^{r,t;f,0) 

xe?(r,0;r',0)gA(r',0;r',t'). 

The fictitious Q will not in general give correct one- 
particle properties. However by definition Q"" gives the 
correct density 

n{r,t) = -2ig<{r,t;r,t) 

(where the factor of 2 comes from spin). Also total cur- 
rents are correctly given by TDDFT. If for instance Jq, 
is the total current from a particular region a we have 



(21) 



Ja{t) = -e I dr —n(T,t) 



where the space integral extends over the region a (e is 
the electron charge). 

The Density Functional Theory and the Runge-Gross 
extension refer specifically to the r basis. However, the 
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arguments remain valid if we instead consider the diago- 
nal density = (clci) in some other basis provided the 
interactions commute with the diagonal density operator. 
The latter condition is essential for the Runge-Gross the- 
orem. Thus, for instance, if the one-electron indices refer 
to a particular lead one can still use Eq. H21|l to calculate 
the corresponding total current (see Section llV)l . 

For later references we now derive an expression for 
the lesser/greater Green function in the linear approxi- 
mation. We consider the partition-free system described 
in the one-particle scheme of Mean Field Theory or 
TDDFT. Let Slf^it) be the small time-dependent effec- 
tive perturbation and SG^'^ = G^'^ ■ SV^ ■ G^^ be the 
first order variation of the retarded and advanced Green 
functions with respect to their equilibrium counterparts 
Go Then, from Eq. ^ we get 



SG^{t;t') 



(22) 



-I- 



dt G^it; t)SU''^{t)G^{0; 0)G^it; 0)G^(0; t') 
dt G^{t; 0)G^(0; f)G^(0; 0)6^^ {t)G^{t; t'), 



where we have taken into account that G^'^ commutes 
with GS=(0;0). The above expression takes an elegant 
form when t' = t. Indeed, for any t > t > one has 
G^{t;0)G^{0;t) = -iG^{t;t) and G^{t;0)G^{0;t) = 
iG^(t;t). Since the integrands in Eq. (|22l) vanish for 
t > t due to the function in G^ in the first term 
and in Gq in the second term, we conclude that for any 
positive time t 

5G^{t-t)^~i j dtG^{t;t)[6U''^{t),G^{0;0)]G^{t;t). 

(23) 

We shall use this equation later on to calculate the linear 
current response in noninteracting and interacting reso- 
nant tunneling systems. 



We also emphasize that all the results of this Section 
remain valid in the interacting case if the bare external 
potential is replaced by the exact effective potential of 
TDDFT, see Section ITTl 

The whole system is described by a quadratic Hamil- 
tonian 



Ho ^ ^ Ekacl^Cka + EqcJcO (24) 

ka 

+ ^ Vka [CfeaCO + clcka] = ^ T^^nci^Cn, 
ka mn 

where a — L, R denotes the left, right lead and m, n are 
collective indices for ka and 0. We assume the system 
in thermodynamic equilibrium at a given inverse tem- 
perature P and chemical potential /i before the time- 
dependent perturbation 

^u(^) = Xl^fe"(^)^i«^fc" + ^o(^)^S'^o 

ka 

= ^ ^ ^m,n(^)Crn^yt 
mn 

is switched on. In principle the time-dependent pertur- 
bation may have off-diagonal matrix elements. In order 
to model a uniform potential deep inside the electrodes 
such off-diagonal terms must be of lower order with re- 
spect to the system size. However their inclusion is trivial 
and it does not lead to any qualitative changes. 

The current from the a contact through the a barrier 
to the central region can be calculated from the time 
evolution of the occupation number operator Na of the 
a contact. From the obvious generalization of Eq. 121(1 
one readily finds 



Ut) = 2e^ 3^[G<fcJt;0 



2eJ2 3ft [G^(t;0)G<(0;0)G'^(0;i) 



O.fco 



(25) 



Vka- 



III. NONINTERACTING RESONANT 
TUNNELING SYSTEMS 

As a first application of the partition-free approach 
we study the time-dependent current response of a non- 
interacting resonant tunneling system. For the sake of 
simplicity the central device will be modeled by a sin- 
gle localized level. All the results of this Section can 
be generalized to the case of a multi-level noninteract- 
ing central device without any conceptual complications. 
There are many different geometries one can conceive 
beyond a one-level model, e.g., a double quantum dots 
model)^ a quantum wire coupled to a quantum dot,^^ 
a one-dimensional quantum-dot arraySi or a mesoscopic 
multi-terminal system^ However, the present paper is 
not intended to give a description of a series of applica- 
tions. Rather, we prefer to illustrate how the partition- 
free approach works in a simple noninteracting model. 



The above expression is manifestly gauge-invariant. In- 
deed, if C/„,„(t) ^ UmAt) + S^n,nX{t) then G"(i;0) ^ 
e-'/ox(t)rftG'^(i.O) while G^{0;t) e'i'o x(*)"'*G^(0; t) 
and the time-dependent shift x{t) has no effect on the 
current response. In the same way it is invariant under 
a simultaneous shift of fi and the initial potential. 
The matrix G^ (0; 0) can be written aa^ 



G<(0;0) 



dC /(C)e 
27rC-5-y' 



where F is the contour surrounding all the Matzubara 
frequencies a;„ = (2n -I- l)m/ (3 + fi clockwise [see Fig. 
in Appendix|2] while rj is an infinitesimally small positive 
constant. It is therefore convenient to define the kernel 



Q„(C;t) - 5] [G^(t;0)G(C)G^(0;i)' 



0,ka 



Vka, (26) 
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with G(C) = [C — £ — V] ^, and to write the current in 
the form 



Jait) = 2e K 



dC 
2tt 



(27) 



It is worth noticing that the partitioned approach 
leads to Eq. ^ with g<(0;0) = if{£) in place of 
G*~(0; 0) = if{£ + V). It is our intention to clarify under 
what circumstances, if any, the long-time behavior of the 
time-dependent current is not affected by this replace- 
ment. 

As a side remark we also observe that Ja{t = 0) in 
Eq. H25|l correctly vanishes. Letting |A„) and A„ be the 
eigenvectors and eigenvalues oiT = £ + V, we have 

J„(0) = -2e^^5[(0|A„)/(A„)(A„|fca)]-t^fe„ =0, 

k n 

since (0|A„) and {Xn\ka) can always be chosen as real 
quantities for systems with time reversal symmetry. 



A. Step-Like Modulation 

The first exactly solvable model we wish to consider is 
a step-like modulation, i.e., Um,n{i) — Q{t)Um,n- From 
Eq. Q it follows that for any t > 

GR(t;0) = -ze~^[^+^+'^l*= / *G^(a;), 

J 27r 

and G^{0;t) = [G^(t;0)]^ The device component of 
G^'^((jj) can be written as 



(28) 



where Eq = eo + Uq. Here, I]«'^(a;) = ^a'^(w) is the 
retarded/advanced self energy induced by back and forth 
virtual hopping processes from the localized level to the 
leads and is given by 



koL 



uj - Ska ± ir] ' 



(29) 



where we have used the short-hand notation Ska — £feo 



1- Uka = Uc,: Steady-State Current 

If Uka = Ua the energy levels of the a lead are 
equally shifted. From Eq. (|25|l it follows that we need 
to estimate the matrix elements G^Q{t; 0), Gq j,^(i; 0) of 
the retarded Green function and the two contractions 
Efe Gl,jO; t)Vk^, Efe G^^, ,J0; t)Vka in the long-time 
limit. We assume that SjJ'^(w) is a smooth function for 



all real u. Then, using the the Riemann-Lebesgue the- 
orem one can prove (see Appendix ^ that the kernel 
QaiC',t) has the following asymptotic behavior 



lim Qa{C;t) = 



de Ta{s) 
de Ta'is) 



G%{e) (30) 



where 

ro(e) 



-23[ER(e)] = 2nY,S{e - Sfc«)Vt. (31) 

k 

In Eq. l|5n|l the r.h.s. has a simple pole structure in 
the C variable and therefore the integration along the 
F contour can be easily performed. Using the identity 
/r(^C/27r)/(C)e''^(C — e)~^ = ifis) the stationary current 



J, 



(S) _ 



hmt 



7(S) _ 



, J a (t) has the following expression 
de rR{e)TL{e) 



where A(e) 



27T [e - eo - A(£)]2 -|- [F(£)/2]2 
x[fie-UL)-f{e-UR)] 
5R[S]" (£)] is the Hilbert transform of F(e) 



jf\(32) 



A{lu) = P 



duj' r{uj') 

271 Ul — to' 



(33) 



It is of interest to note that the dependence on the bias 
Ua appears not only in the distribution function / but 
also in the quantities F and A, see Eqs. (|31|I - H33I) . The 
dependence of the self energy on the level-shifts is phys- 
ical since when the particle visits the reservoirs experi- 
ences the applied potential. We also remark that Eq. 
(|32|l is of the Landauer typef2£ More generally the Lan- 
dauer formula is valid for any mesoscopic device pro- 
vided it is noninteracting. This result agrees with the 
one obtained in the partitioned approach by Jauho and 
coworkers There the leads are decoupled from the 
central device and in thermal equilibrium at different 
chemical potential fi^ and fj,R, and inverse temperature 
(3l and Pr in the remote past. In order to preserve 
charge neutrality each energy level eka must be shifted 
by /ia — /i where fj, is the chemical potential of the two 
undisturbed leads. The stationary current is then ob- 
tained by switching on the contacts, i.e., the hybridiza- 
tion part of the Hamiltonian. By tuning Pr = (3^ — P 
and hr — fJ-L — Ur — Ul the current is given by Eq. 

To summarize we have found that for noninteracting 
leads a steady state develops in the long-time limit when- 
ever 1) The one-body levels of the charge reservoirs form 
a continuum and 2) The self energy due to the hop- 
ping term is a smooth function. Under these hypothe- 
ses the time-translational invariance is restored by means 
of a dephasing mechanism. The comparison of our re- 
sult with the one obtained in the partitioned scheme 
provides the criteria of equivalence: besides the tuning 
f^L — = Ul — Ur one needs to shift the levels of the 
a reservoir by fia — fJ-. 
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2. Uka ~ Uc,: Time- Dependent Current in the Wide-Band 

Limit 

The calculation of the stationary current is greatly sim- 
plified by the long-time behavior of the various terms 
coming from Eq. I|25|) . However, as far as we are inter- 
ested in the current at any finite time we need to spec- 
ify the structure of the retarded (advanced) self energy. 
Here, we consider the so called wide-band limit where 
the level-width functions Ta{Lo) = 27q, are assumed to be 
constant and hence, from Eq. (j32J), Ka{Lo) = 0. In this 
case G^q{u)) has a simple-pole structure and the calcu- 
lations are slightly simplified. We emphasize that what 
follows is the first explicit result of a time-dependent cur- 
rent in a model system in the framework of a partition- 
free approach and therefore also a simple model could 
be of some interest. Without loss of generality we can 
always choose Eq = 0; for the sake of simplicity we also 
consider Uq = Q. We defer the reader to Appendix O for 
the details. Here, we just write down the final result for 
Ja(i): 

J„(i) = -467.6-^*1 ^/H (34) 



C/a'e"^* + 2tJC0s[(w + Ua')t] + 27sin[(tj + Ua')t] 



where Jq^^ is the stationary current of Eq. I|32|l and 
J = IR + JL- One can easily check that 1) For t — > oo Eq. 
(EH) yields the result in Eq. 2) For t = the cur- 

rent vanishes, that is Ja{0) — and 3) For Ul — Ur — 
the current vanishes for any t. Eq. H34|l can be rewrit- 
ten in a more physical and compact way if we exploit 
the particle-number conservation. Denoting by ng the 
particle number operator in the central device we have 



which are absent if the system is uncontacted for nega- 
tive times. We have explicitly verified that by discarding 
these terms our formula reduces to the one obtained in 
the partitioned scheme. For long times, the extra terms 
vanish and our scheme reproduces the earlier steady-state 
results. 

If one of the two leads does not undergo any level shift, 
e.g. Ur = 0, from Eq. we get 



'It 



duj 



C/ie"^* -h 2wcos[(w -I- UL)t\ -h 27sin[(a' -f- C/l)^] 



(35) 



The transient behavior of the time-dependent quantity 
Ja{i') — Jq is not simply an exponential decay. In Fig. 
121 we have plotted JR(t) in Eq. (|35|l versus t for differ- 
ent values of the applied bias at zero temperature. 
The current strongly depends on Vh for small Vh while 
it is fairly independent of it in the strong bias regime; 
using the parameter specified in the caption, the time- 
dependent current has essentially the same shape for any 




FIG. 2: Time-dependent current Jn(t) for different values of the 
applied bias Uj^ = 0.8, 2.0, 4.0 and 6.0. The numerical integration 
has been done with 'fR = fL = 0.2, fi = and zero temperature. 



so that 



7 ^^-\ 7(3) , d 

Mt) = Jr +^—-^(^"' 



4e2^e-^* 



duj 
2^ 



Ur 



UJ + Ur + i'y 



-Ul 



UJ + Ul + i'y 



JL{t) is obtained by exchanging i? <~> L in the r.h.s. of 
the above expression. Therefore, JR{t) ^ ~JL{t) for any 
finite time t, even in the symmetric case 7_r = 7l; the 
time derivative of (no) contributes to Jr and Jl in the 
same way. 

Our formula for the nonlinear transient current clearly 
differ from the one obtained by Jauho et al& in the parti- 
tioned scheme. Indeed, the prescribed integration along 
the imaginary axis gives extra terms (see Appendix Q 



7=2.0 




y= 1.0 



7=0.4 



1 2 3 4 5 6 7 



FIG. 3: Time-dependent current for Ul = 6, = and zero 
temperature versus time for three different values of the line widths 
7^ = 7^= 0.2, 0.5 and 1.0. 



In Fig. 121 the current Jfl(t) is plotted for different val- 
ues of the total line width 7 and for a fixed value Ul = 6 
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of the applied bias. As expected, the larger is 7 the bigger 
is the slope of the current in t = 0. 



0.5 y=2.0 




FIG. 4: Stationary current versus the applied bias at zero temper- 
ature and chemical potential for three different values of the line 
widths 7i? = 7l = 0.2, 0.5 and 1.0. 

Finally, in Fig. ^we report the trend of the stationary 
current J)j as a function of the bias Ul for three different 
choices of the level widths. As one can see the bigger is 
7 and the wider is the range of validity of the Ohm law. 

B. Arbitrary Modulation: Theorem of Equivalence 

We have shown that the steady-state current induced 
by a step-like modulation does not change if one uses the 
of the partitioned approach, given by the first term 
on the r.h.s. of Eq. H14|) . in place of the one coming 
from the partition-free approach. This reasonable result 
is now proved and not simply postulated. The equiv- 
alence between the two expressions for the current is of 
special importance since it is much easier to work in the 
partitioned scheme. However, it has been proved only 
for step-like modulations with Uka — Ua- Here, we prove 
that the above equivalence remains true under very gen- 
eral assumptions. To this end we consider the quantity 

s^;v(^;i') = E^r<f(^;i')vL (36) 

k 

where V is an arbitrary complex function of ka. Then, 
the following theorem holds 
Theorem of Equivalence: If 

lim J:%{t; t') = hm (i'; t) = (37) 

for any nonsingular V, then 

lim[Q„(C;0-?a(C;i)] =0, (38) 

t—*oo 

where qc{C,t) = Efcl^^"^ 0)g(C)G^(0; t)]o,fcal4a, and 
s{C) = [C ~ S]~^ is the uncontacted Green function. 

Eq. (|38|) says that if we apply the same time-dependent 
perturbation the same asymptotic current will emerge in 
the partitioned and partition-free approaches. 



Proof: In terms of the self energy Ey = ^aV^ 
equation of motion for G^q takes the form 

- eo(o} G^oif. t') - [^v ■ GU{t; t') = 5{t - t% 

where the symbol • denotes the real-time convolution. 
We now consider the limit t ^ 00. The hypothesis H37|) 
implies 

lim G^o(t;t') = 0, (39) 

t — *oo ' 

which in turn implies that [Ey • G^f^{t]t') 0. Fur- 
thermore, from the Dyson equation for G^^.^, wc find 

Y.<^o.{t'-.t)Vk^ = [G^,o ■ <v](^';t) 0. (40) 

k 

We note that the above two asymptotic relations have 
been obtained for i' = in the special case of a step- 
like modulation, see Eq. ljBl|) . Here, we have shown that 
they hold in a more general context. As a consequence of 
these two identities, the asymptotic difference [QaiC'j ^) ~ 
^(Cii)] can be written as 

[Glo ■ (t; 0) X Go,o(C) x ■ [6 + G^q • J^^v] (0; t). 

Here, E^} = J^a ^a.v is given by Eq. ^ with = 
Vka/iC - £fea). Since C G T, 3[C] ^ and hence V^^ 
is nonsingular, meaning that Eq. H37I) holds. Eq. I|37() 
together with Eq. 14U|I imply the equation of equivalence 

As a simple application of the Theorem of Equivalence 
one can calculate the stationary current for an arbitrary 
step-like modulation. The quantity qa{C]t) is simply 
given by the first two terms of Eq. (|B2p . Both have a sim- 
ple pole structure in the ( variable and we can perform 
the integration along the contour F. Using the definition 
in Eq. H31|) . with i^a — £ka + Ua{ska), one obtains 

4' - -ej g/(e) (41) 
X {Ff (£)F«(e + ULie)) |G^„(e + UUe))]' 
+Tfis)TUs + Un{s)) \Gl,is + Unis))\'} 

The quantity rL°^(e) ee 2T:J2k ^ka^i^ ~ ^ka) is the equi- 
librium line width. Eq. 141|l reduces to Eq. (|32|) if 

Ua{£) — Ua since in this case Ta\e — Ua) = Ta{£)- 

In the noninteracting case it is reasonable to assume 
that Eq. H41I) yields the steady-state current even 
for an arbitrary time-dependent disturbance such that 
\\mt->ooUka{t) = Uka and limj^oo C/o(i) = Uq. In the 
next Section we shall prove that the asymptotic current 
has no memory and depends only on the asymptotic value 
of the external perturbation. 
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C. Memory Loss Theorem 

If the condition (|37|l of the Theorem of Equivalence 
is fulfilled, the asymptotic value of the nonlinear time- 
dependent current in Eq. H27|) simplifies to 



■/J/(C)e''^..(C; 



t) 



(42) 



2e 3? 



kl3 



k' 



We note in passing that expressing Gq ^.^ and G^^ j,,^ 
in terms of Gq g and G^g respectively, Eq. can be 

V e-< 1/2 

oka ^ k(x 



rewritten in terms of y 



Ut) = 2e 5ft{[G;^g • S<y](i;0 + [Gg<g • ^ivW^t)} , 

where the asymptotic relation Ggg — ^o,o ''^aV ' 
G^g has been used [see Eq. (j^OI]- This agrees with the 
result obtained by Wingreen et alS^ in the partitioned 
approach, as it should. 

In general Ja{t — > oo) is not a constant unless the 
external perturbation tends to a constant in the distant 
future. In this case the following theorem holds 
Memory Loss Theorem: If 

lim Ua{e,t) = C/a(e), lim [/o(t) = Uq 

t — *CXD t — >OC 

the current Ja(t) tends to a constant, given by Eq. 
in the long-time limit. 

Proof: Is convenient to denote with G and g the Green 
functions corresponding to the step-like modulation with 
coefficients Ua{£) and Uq. We have already shown that 
in the long-time limit Eq. (05)1 yields Eq. ||tl} if G^'^ = 
(JR,A^ The Memory Loss Theorem is then proved if 



iim — „ , 



lim 



Y.k' Gtl3,k'aiO'^tWk'a 

Ek' Gtf^,k'Mt)Vk'a 

(43) 

for some real constant A^^ . 

According to Eq. (|^ . the device component of 
the retarded Green function Gg'g(t oo;t') vanishes 
for any finite t' . Since \imt,t'->oo[go'it;t')/g^{t;t')] = 
1, from G|^g = g^ + go ■ ' '^o;o follows that 
limt^t'^oo[GQ'Q{t;t')/G^Q{t;t')] — 1. Let us now consider 
the Dyson equation 



Glkcit;t') = / diGloit;mo.gUt;t'), 



(44) 



with t' = and t ^ oo. Since the integrand vanishes for 
any finite t, we can substitute Gg^g with Gg^g. Further- 
more, since the applied bias tends to a constant in the 



distant future, \imt^^[gf^{t;0)/gf^{t;0)] = e-*^'^" for 
some real quantity Aka- The l.h.s. of Eq. H43|) is then 
proved. A similar reasoning leads to the r.h.s. of Eq. 

63- 



D. Linear Response in the Wide-Band Limit 

In the case of small time-dependent perturbations, one 
can use Eq. H23|l to calculate the lesser Green func- 
tion. In order to carry on the calculations analyti- 
cally we consider the wide-band limit and we choose 
5Uka{i) — 5Ua{t). For simplicity we omit the subscript 
in the retarded and advanced equilibrium Green func- 
tions. By explicitly writing down the matrix product in 
Eq. H23|) one readily realize that we have to calculate the 
functions G«:g(t - t),Gl^,^,{t - t), Efc VkaG^^,Jt-t) 
and VkaG^,^, kai^ — t). They are easily obtained from 
Eqs. HG2|) by simply replacing e — > e and t ^ t — t. 
The calculations are rather similar to those already per- 
formed to derive the expression H34|) and they are left 
to the reader. Denoting by SJa{t) the time-dependent 
current in the linear regime one ends up with 



SJa{t) = 4e7„sR J / dt 



duj gi(w-Oo)(t-t) 

T^fi^) a 



5Uo{t)-SUa,{t) + 2iJ2lc.^ 



SUo{t)-SU^,{t) 

ClJ — ilg 



(45) 



where r^o — eo — ij. In the special case Eq — 0, SUo{t) = 
and 5Ua{t) = SUa = const, 5Ja{t) reduces to the time- 
dependent current in Eq. 1)34(1 to first order in 5Ua , as it 
should. 

Eq. takes a very simple form = 5Uq = 5Un = 
and a — R: 



5JR{t) = 4e7fl7L / dtSUL{t-t) 



S[/(n)] „-2.f 



7 



f3 



E 

.n=0 



g(ii^„-7)t 



where LUn — {2n+ l)TTi/P + fj, are the Matzubara frequen- 
cies and the identity 



doj 

2^^(") 



^2 _[_ j2 



-y 



27 /3,t^g-^+7^ 



has been used. In the special case of a vanishing chem- 
ical potential, the Matzubara frequencies are imaginary 
numbers and 6Jji{t) simplifies 



6JR{t) = -2e 



7fl7L 
7 



dt6UL{t~t)e-^^* 



x|tan[^] 



-[f-7]f 



e ^f* 



-m 



(46) 
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where 



Fit) - E (-)"'i>[e-^,l 



i ^ 7r+(-r/37. 



m=0,l 



27r 



is a linear combination of the Lerch transcendent func- 
tions $[z, s, a] = X^^^o ^"/(a + "-)"■ 



\ p=ioo 

S=2 



(a) 




P=700 



(3=2 



(b) 



FIG. 5: Current versus time for two different external disturbance 
SUl- In both figures ■jr = 7l = 0.5; the current is plotted for two 
different inverse temperature /3 = 2 and /3 = 100. In (a) SU^ is a 
square bump-like modulation whose duration is 1 while in (b) the 
duration is 5. 



In Fig. 13 we show the trend of SJii{t) for square bump- 
hke modulations. On the top SUL{t) = 9(^)6(1 - t) 
while on the bottom 5UL{t) — Q{t)Q{5 — t); both distur- 
bances are considered for two different inverse tempera- 
ture (3 = 2 and /3 = 100. As one can see the effect of 
an increasing temperature consists in a sort of rescaling 
of the time-dependent current. The line widths 7^ have 
been taken equal and large enough to justify the linear 
approximation. Since the disturbance is of order 1, from 
Fig. ^one can see that 7_r = 7l = 0.5 is a good choice. 

The ac current in the linear approximation is plotted 
in Fig. for /3 = 100 and 7^ = 7l = 0.5. The time- 
dependent disturbance is taken to be 5UL{i) = sinwo^ 
with Wo = 5, 10, 20 in (a), (b) and (c) respectively. 
Finally, in Fig. E|we have considered the current response 
to a periodic square bump-like modulation for different 
values of the period. 



57„ 

0.04 
0.02 
0.00_ 



o)„=5 



(a) 



o)„= 10 




FIG. 6: Current versus time for an oscillating external disturbance 
SUl- In both figures = 7l = 0.5 and the inverse temperature 
is /3 = 100. SUL{t) = sin wo* with ojq = 5, 10, 20 in (a), (b) and 
(c) respectively. 




IV. INTERACTING SYSTEMS 

In earlier theoretical works on quantum transport one 
can distinguish at least two schools. In one school one 
tries to keep the full atomistic structure of the conductor 
and the leads, but all works so far are at the level of the 
Local Density Approximation (LDA) and only the steady 
state has been considered. The advantage of this ap- 
proach is that the interaction in the leads and in the con- 
ductor are treated on the same footing via self-consistent 
calculations on the current-carrying system. It also al- 
lows for detailed studies of how the contacts influence the 
conductance properties. 



FIG. 7: Current versus time for three different periodic square 
bump-like modulations. In any figure ^ji = 7l = 0.5 and j3 = 100. 
The thin lines represent SUl while the thick lines represent SJ^. 



The other school is using simplified models which al- 
lows the analysis to be carried much further. Consider- 
able progresses have been made in this respect for a lo- 
calized level described by a Lundquist-like model^ii^SiS 
and for the so called "Coulomb island"i^^i^ where Hq in 
Eq. (|24|) is replaced by the Anderson Hamiltonian. How- 
ever, all these works treat the leads as non-interacting, 
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which prohibits a reahstic description of the contacts and 
of the long-range aspects of the Coulomb interaction^. 
The model approach are based on a partitioned scheme 
which makes the time-dependent results difficult to in- 
terpret. 

We here want to show how the current LDA by Lang 
et al}^'^^ follows from the TDDFT scheme described in 
Section 111 CI We also present an exact result for the 
steady-state current of an interacting resonant tunneling 
system. Finally, the transient behavior of a capacitor- 
device-capacitor system is investigated on the level of 
Mean Field. 



A. Steady-State Limit of TDDFT 

In Section lnil we showed that under certain conditions 
a steady state is reached in the long-time limit, and that 
this limit is independent of history. We also showed 
that the partitioned and partition-free treatments give 
an equivalent description of the steady state. The mech- 
anism for the loss of memory was pure dephasing, and it 
holds provided the leads are macroscopic while the device 
is finite. Another important ingredient is that the applied 
bias is uniform deep inside the leads. With these assump- 
tions, our results can be generalized also to more general 
cases than the simplified model explicitly considered in 
Section |llll In TDDFT, the full interacting problem is 
reduced to a fictitious noninteracting one and all the re- 
sults of Section imi can be recycled. In the case of Time 
Dependent Local Density Approximation (TDLDA) , the 
exchange-correlation potential Wxc depends only on the 
instantaneous local density and has no memory at all. If 
the density tends to a constant, so does the effective po- 
tential If^^ , which again implies that the density tends 
to a constant. Owing to the non-linearity of the problem 
there might still be more than one steady-state solution 
or none at all. 

If a steady state is reached in TDDFT, we can go di- 
rectly to the long-time limit of the Dyson equation and 
work in the frequency space. We may with no restriction 
use a partitioned approach and split the fictitious one- 
electron Hamiltonian matrix in a non-conducting part £ 
and a correction V involving one-body hopping terms 
between the two leads and the device. The lesser Green 
function of TDDFT fulfills 

e<(£) = [1 + g<(e)[l + Vg^{e)], 

where g is the uncontacted TDDFT Green function [cf. 
Eq. (|13|l ]. In direct space, the imcontactcd g< can be 
written 

5<(r, r',£) 2ni^ fa{ema)<t)7na{r)(j)*„iail^')S{e - Cma) 
ma 

in terms of diagonalizing orbitals (fima with fictitious 
eigenvalues Cma for the left and right leads (a — L,R) 
and the device (a = D) and Fermi functions fa with 



chemical potential /Iq. The chemical potentials for the 
two leads differ, and the final result is independent of the 
chosen chemical potential for the device. When we apply 
1 -t- g^V = g^[g^]~^ to an unperturbed orbital (^„ia, 
it is transformed to an interacting, i.e., contacted eigen- 
state ipma- Above the conductance threshold, states orig- 
inating from the left lead become right-going scattering 
states, and states from the right lead become left-going 
scattering states. In addition, fully reflected waves and 
discrete state may arise which contribute to the density 
but not to the current. Thus, 

C/<(r,r',e) = 27rz ^ /a(e™a)V'ma(r)V'ma(r')'5(£ - 

These results correspond closely to the general approach 
by Lang and coworkers i^^ii^ In their approach, the con- 
tinuum is split into left and right-going parts, which are 
populated according to two different chemical potentials. 
The density is then calculated self-consistently. Lang et 
al. further approximate exchange and correlation by the 
LDA and the leads by homogeneous jellia, but apart from 
these approximations it is clear that his method imple- 
ments TDDFT, as described in Section lll CI in the steady 
state. It is also clear that the correctness of Lang's ap- 
proach relies on the Theorem of Equivalence between the 
partitioned and partition-free approaches and the Mem- 
ory Loss Theorem derived here. The equivalence between 
the scattering state approach by Lang et al. and the 
partitioned non-equilibrium approach used by Taylor et 
al.i^ii^ has also been shown by Brandbyge et al.^"^ 

As shown above, the steady state of TDDFT can al- 
ways be formulated in terms of orbitals which diagonal- 
ize the asymptotic one-particle Hamiltonian matrix. The 
current-carrying orbitals can always be grouped into a 
right-going class and a left-going class. As a consequence, 
the current can be expressed in a Landauer formula 

Jr'' = -e'YyfL{emL)TraL " fRiemR)%nR] (47) 
m 

in terms of fictitious transmission coefficients 7^„a and 
energy eigenvalues ema, a — L,R. We also wish to em- 
phasize that the steady-state current in Eq. H47() comes 
out from a pure dephasing mechanism in the fictitious 
noninteracting problem. The memory-loss effects from 
scatterings is described by ^xc and fxc- 

B. One-Level Resonant Tunneling System 

In this Section we consider a resonant tunneling system 
described by the quadratic Hamiltonian of Eq. H24|) and 
an inter-particle interaction 

Hw — 2 ^ Wm^nnmnn, 

where = c^^Cm is the occupation number operator of 
the level m and Wm,n = Wn,m, is a symmetric matrix. (If 
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Hw includes long-range terms, the regrouping of poten- 
tial terms as discussed in Section Hi CI must be done.) In 
the generalized TDDFT scheme (based on the rim occu- 
pations rather than on density) outlined in Section III CI 
the fictitious Green function Gm.n is obtained by solving 
the Dyson equations with K = T + U'^^ , where 

it)]. 

If t7°^ satisfies the hypothesis of the Theorem of Equiv- 
alence and of the Memory Loss Theorem we can use Eq. 
(|41|l and write an exact formula for the steady-state cur- 
rent of an interacting resonant tunneling system: 



(S) 



de 
2^ 



fie) 



(48) 



X {r^°\s)rjiis + Uf{e)) \gl,{e -f Uf{e))\^ 

+Tf{e)TLie + Uf{e)) \gl,{e + <(e))f } . 

For normal-metal electrodes we expect that the effec- 
tive potential U^{e,t) U^^{e) — const provided 
Uaie,t) Uaie) — const when i — cx). The con- 
stant Uf{e) may depend on the history of Uais,t) 
while the steady-state current is independent of the 
history of Uf{e,t). Q^^^iuj) is given by Eq.||2Hl) with 
£o = £o + linit->oo U^^{t) and with Yp" from Eq.lj^HIl with 
£ka — £ka + hmt_>oo U^^{t). For the sake of clarity, Eq. 
(I48|l has been written for systems having a one-to-one 
correspondence between the one-body indices ka and the 
one-body energies eta + U^ai^)- ^he generalization to 
systems with degenerate levels is straightforward and it 
is left to the reader. 

As a further example we study the RPA time- 
dependent current response in the partition-free ap- 
proach. In the Hartree approximation the Green function 
satisfies the equation of motion (|18l) with Sc = and 

where nf-iz) — — iG'^;^(z; z). According with the results 

obtained in Sectional the lesser Green function G^'^ 
is given by Eq. ^ with G ^ G^. Therefore, in the 
linear approximation we have 

6G<{t;t) (49) 

= -iJ dtGo'^{t;t)[SU''^{t),G^-<{0;0)]Go'^{i;t), 



with 



sw^{t) = su{t) + si:^{t). 



(50) 



Eqs. H49|I - H5U|I form a coupled system of integral equa- 
tions for the unknowns dG^{t;t) and 6U^^{t). For a 
capacitor-device-capacitor system one can take 







if m = ka, n 
otherwise 



k'a' 



Thus, putting an extra particle in the isolated a capacitor 
costs an energy Waa per particle. This means that the 
transfer of a finite number of particles from one capacitor 
to the other causes a finite change of the effective applied 
bias. We expect that the current vanishes in the long- 
time limit unless the applied bias continues to grow up. 
The coefficients Wrl = Wlr mimic the repulsion energy 
between two particles in different capacitors. Actually, 
one can also consider the interaction between a particle 
in the central device and another in one of the two ca- 
pacitors. No extra comphcations arise if Wo.fca = Wa, 
Vfc, and the results we are going to obtain can be easily 
extended. 

Switching a bias SUkait) = SUait), from Eq. (|50|l one 

gets SU^Jt) - Sm^JUfit) with SU',^it) - SUfit), 
Vfc, and 



5Uf{t) = 5U^{t) 



i V / diWo.f,5Mt), 

^ a Jo 



(51) 



where it has been taken into account that SN^(t) = 

T,k^^fait) = --efodt SJc^it). Since SW^ has the 
same matrix structure of the bare SU, in the wide-band 
limit the linear time-dependent current 5 J ait) is given 
by Eq. with 5U replaced by . (It is worth 

noticing that the wide band limit still makes sense if the 
line width is approximately constant in a small interval 
around the chemical potential fi.) In this way the sys- 
tem of Eqs. (|l^ - lt5n|l is reduced to a system of 4 coupled 
integral equations for the 4 scalar unknowns SU^^, SJ^ 
with a — L,R. The symmetric case — = 7/2, 
Wrr = Wll allows a further simplification. Let us de- 
fine 5Uf = 6Uf ± SUf, 6U± = 6Ur ± 5Ul, SJ± = 
5Jr±5Jl and W± = Wrr±Wrl = Wll±Wlr. Then, 
from Eq. (|51|) we find 



5Uf{t)^5U±{t)- — 



[ dtdJ±{t), (52) 

^0 



while from Eq. 1)45(1 



dJ+{t) = 2e7 / dt C+{t - t)[25Uo{t) - 6Uf{t)], (53) 
Jo 



SJ-{t) ^ -2e7 / dt C-{t - t)SUl^{t), (54) 
Jo 

where 



C±it) = ^i I pfiu;)'- — 



is the conductivity kernel. Once 6J± has been obtained, 
one can calculate 5Jr = {5J++5J-)/'^ and &Jl = (<5J+ — 
(5J_)/2. 

In order to illustrate what is the time-dependent re- 
sponse of this model we have considered the zero temper- 
ature case with SUR{t) = -SULit) = (l/2)e(t)e(io - t) 
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5/ 





FIG. 8: Numerical solutions of Eqs. l52i - iHi^ - l54i in the zero 
temperature limit with fi = eo = 0, W— = 5 and an external 
disturbance as described in the main text. The thick lines are the 
current in (a) and the effective potential in (b) for the step like 
modulation. The currents and the effective potentials for to = 
1, 2, 4, 6 unstick from the thick line and start to oscillate and 
eventually vanish after a time t oc 1/W—. The vertical lines are the 
bare applied potentials. 



and 5UQ{t) = 0. Then, 5U+{t) = and hence 5J+{t) = 
SUf{t) = 0. It foUows that SJait) = -SJL{t) = 
SJ-{t)/2 and 5Uf{t) = -SUf{t) = 6Ul^{t)/2 for any 
time t. In Fig. |Sfa) we display the time-dependent 
current for square bump-hke modulations with = 
1, 2, 4, 6 and W- ~ 5. The thick line is the current 
for the step like modulation SU-{t) — 8(t); depending 
on the value of to the current unsticks itself from the thick 
line giving rise to different damped oscillating curves. In 
correspondence of each to a vertical line has been drawn; 
it represents the bare applied potential 6U-{t). Fig. |HJb) 
shows the time-dependent effective potential 6U^{t). As 
the current response, it drops to zero in the long-time 
limit since the interactions completely screen the applied 
bias after a time t oc 1/W-. 



V. SUMMARY AND CONCLUDING REMARKS 

In the present work we have used a partition-free 
scheme in order to treat the time-dependent current re- 
sponse of a mesoscopic system coupled to macroscopic 
leads. To this end, we have further developed the 



Keldysh formalism and we have formulated a formally ex- 
act theory which is more akin to the way the experiments 
are carried out. Among the advantages of the partition- 
free scheme we stress the possibility to calculate physical 
dynamical responses and to include the interactions be- 
tween the leads and between the leads and the device in 
a quite natural way. 

In the noninteracting case we have shown that a perfect 
destructive interference takes place provided the energy 
levels of the leads form a continuum. The steady-state 
develops due to a dephasing mechanism. The comparison 
of our steady-state current with that obtained in the par- 
titioned scheme shows that the two currents are equiv- 
alent if the energy levels are properly shifted in order 
to preserve charge neutrality. This kind of equivalence 
remains true for any time-dependent external potentials 
(Theorem of Equivalence). The Theorem of Equivalence 
has then been used in order to prove that the steady- 
state current depends only on the asymptotic value of 
the external perturbation (Memory Loss Theorem). For 
the sake of clarity, the Theorem of Equivalence and the 
Memory Loss Theorem have been proved for a single-level 
central device. The generalization to a multi-level central 
device is straightforward, as can be readily verified. In 
the wide band limit we have obtained an analytic result 
for the time-dependent current in the case of a step-like 
modulation and for arbitrary modulations in the linear 
regime. 

The interacting case represents a more difficult chal- 
lenge and the expression for the lesser Green function 
at any finite time is more complicated than that com- 
monly used to calculate steady-state response functions. 
As an alternative to a full many-body treatment we have 
proposed a formally exact one-particle scheme based on 
TDDFT. Then, all the results obtained in the nonin- 
teracting case can be recycled provided we substitute 
the external potential with the exact effective potential 
of TDDFT. Although it is difhcult to prove any rigor- 
ous results for the effective TDDFT potential, we expect 
the interactions to reduce the memory effects even fur- 
ther compared to the noninteracting case. Thus, any 
nonlinear steady-state current can been expressed in a 
Landauer-like formula in terms of fictitious transmis- 
sion coefficients and one-particle energy eigenvalues. The 
steady-state current depends on history only through the 
asymptotic shape of the effective TDDFT potential. This 
exact result may prompt for new approximations to the 
exchange-correlation action functional A^^c- In the effec- 
tive one-particle scheme of TDDFT the steady-state cur- 
rent comes out from a pure dephasing mechanism. The 
damping mechanism (due to the electron-electron scat- 
terings) of the real problem is described by A^c- As an il- 
lustrative example we have also calculated the RPA time- 
dependent current of a capacitor-device-capacitor system 
and we have displayed the effect of the charge oscillations 
in the discharge process. 
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APPENDIX A: PROOF OF EQ. llToll 

It is convenient to define Go as the solution of Eqs. 
(|18|1 with Sc = 0- Gq satisfies all the relations we have 
derived for a noninteracting system in the presence of an 
external disturbance. By using the Langreth theorem, 
we get 



<5 + G« ■ 



G^-S^ + Gl*sf 



+G^ ■ Sl ★ G|i + gU S ★ g|^ 

and solving for G^ 

G^ = [(5 + G^ • S^] • Gf • [5 + • G^;^] 
+ ■'E^ + G^ ★sH -G^ 

+[G" • Sl * g|^ + gU £ ★ g|^] • [(5 + • G^]. 

Next, we use 

Gf(t;i') = G|^(t;0)Gf(0;0)G^(0;O 

and 

Go(r;i) = -iGo(T;0)G^(0;i), 

so that 



G^{t;t') = G^(i;0)Gf (0;0)G^(0;t') 



G« . • G^ 



gUs^.g' 



(Al) 
{t;t') 



G^ • Sl * Go + gU S ★ Goj {t; 0)G^(0; t'). 

As in the noninteracting case, we proceed by writing 
down the Dyson equation for G(t;T). Taking into ac- 
count that 

G(r; r') - Go (r; r') + [Go * S * G] (r; r') (A2) 
and that 

Go(t;r)=zG?(t;0)Go(0;r), 

we have 

G{t; t) = \g^ ■ * gI {t; r) + iG^{t; 0)G(0; r). (A3) 



Similarly, it is straightforward to show that 

G(T;t)= [g^sT-G^I (T;t) -iG(T;Q)G^{0;t). (A4) 

Substituting Eq. (|X3l) into Eq. IjAip and using Eq. ((X2l) 
one finds 



G'^{t;t') ^ G^(i;0)G^(0;0)G^(0;t') 
+ \g^ • [s?^ + s1 *G*Sr 



G' 



{t-t') 



iG^{t;0) [G*Sr.G^j {0;t') 

G^-Sl *Gj (t;0)G^(0;t')- (A5) 

Using Eqs. (|A3|I - (|A4(I to express the last two terms as 

iG^{t;0) [G*Sr.G^j {Q;t') 
= ^G"(^;0)G>(0;^') -G^(i;0)G>(0;0)G^(0;t') 

and 



G^'' ■ E I ★ G 



{t;0)G^{0;t') 



= -iG<{t;0)G^{0:t') - G^{t;0)G<{0;0)G^{0;t'), 
we end up with Eq. (|19|l . 

APPENDIX B: PROOF OF EQ. I^Ot 

Due to the smoothness of the self energy, in the long- 
time limit we can use the Riemann-Lebesgue theorem to 
obtain the following asymptotic behaviors 



hm G;^o(^;0) - lim ^G^,j0;t)l4„ = 

k 



(Bl) 



and 



hm Gl,^{t;0) = -*yfcae---*Gj^o(efcJ. 



\imY,Gta',ka(.0;mc.=iVk'c.'e''^'-'' 

k 

X [<5a,a' + GQQ{ik'a')^ti^k'a')] ■ 

From the above results and the definition (|26|l one has 
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lim Q4C;t) = ^ -J^G«_„(£fe,„) 



(B2) 



v3 , 



Z2 /■ C!o'.oi^k'a')GQQ{ek'a')^ai^k'a') 
k'a' 



lim Go,o(C) E 



7i C - £fe'. 



- lim Go,o(C) E G^o(^~fe'o')e-'^~'='"'* 



1/2 

J.// C ~ £fc"Q 



T/2 



Gq o(£fc"ct")^a(^fc"a!")6 



r„(£) 

+ lim Go,o(C) / ^T^{e) 



+ ^lim Go,o(C)E y i^r„,(£')G^o(£' 

a' 

/rip" 
^-V^„{e")GUe")^t{e") 



de' . . e-'^'* 



C - £' + f/a 

C - e" + C/a" 



where the relation 



^ Ska.k'a' Vka ^k' a' ^ , 
'~^ka,k'a'[(,) — . _ ^ H "7^ "T^ (-fO,ol<,j, 

has been explicitly used. Since C £ T the quantity [C — 
£ — Uci]~^ is a smooth function of e for any real e. In 
the limit t ^ oo the last two terms in Eq. ljB2p vanish 
according with the Riemann-Lebesgue theorem and Eq. 
(|5n|l is recovered. 



APPENDIX C: PROOF OF EQ. 

The quantity Qa (C; t) involves the multiplication of 
three matrices and we can recognize four contributions, 
two containing Gq g and other two containing Gq . It 
is straightforward to verify that 



1 o 



3[C] > 



Go,o(C) 



(Cl) 



and that G^'(^{uj) = [lu ± ij] ^, where 7 = 7_r + 7l- 
Hence 



G?:o(^;0) = -^e-^^ 

£ka + «7 

E^fe"^^fc"(0;i) = -7ae-^*, 
fc 

Y.^k'n'Gt.,,'o.'{0;t) = iS^.,,,Vkae" 

k' 

-la'Vko 



(C2) 



Ska - »7 



C plane 

r=r+u r_ 



3iti/p 

-it</p 

-3it //p 
-5it «7p 



FIG. 9: Contour T of Eq. j^. The black dots correspond to the 
position of the Matzubara frequencies in the complex C, plane. 



Eqs. (|C1I) - (|C2|) are all what we need in order to 
evaluate the quantity Qa{C'-,t) in Eq. (|26() . The 
time-dependent current is then obtained integrating 
QaiC', t)f{C)^^^ over C along the contour F of Fig. |51 ac- 
cording with Eq. H27|l . Using Eqs. ljC2p and expressing 

GfcQ,o(C) = Go,fcQ(C) and GkaM'a'iC) in terms of Go,o(C) 
we obtain 

G«,o(i; 0)Go,o(C) E VkcGl.^iO; t) = z7aGo,o(C)e-'^*, 

k 

(C3) 

J2G%{t;0)Go,k'a.'{C)J2^k»Gt^,,U0;t) = (C4) 



= 27„e-^*Go,o(C) 
2i7ae-^*Go.o(C)ET" 



de 



27r C - e + C/c 
de e*=* - e-'^* 



-I- 



27r (C-e + C/„0(e-«7) 



We are left with the contributions containing G^^,^,. 
One of them is quite easy to evaluate and yields: 

EG«:,,„,(i;0)Gfe.„,,o(C)E^fc"G^fc"(0;t) = (C5) 



= 2i7aGo,o(C)E'^«' 



de e"*^*-'^* - e-^"^* 
2^(C-e + C/„0(e + «7)' 
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The other one is much more involved, but nothing more 
than standard algebra is needed to get the following ex- 
pression 



E E ^J^fc'"' 0)Gfc'„',fc"a" (C) X (C6) 



k'a' k" a." 



de 



1 - e^^'e""*"* 



27r (e + i7)(C-£ + C/a) 



+ 



2«7a E ^ 
a' 

2zGo,o(C)E^"' 



1 



27r C - e + C/a' 



e + «7 



de e"*^* - e-'^* f de' 
2^(e + i7)(C-e + C/c.O i 2^ 



^^ie't _ p — ft 



C-e' + C/„ (e'-i7)(C-£' + C^a") 



The r.h.s. of the above four equations must now be multi- 
plied by /(C)e''^ and integrated over C, along the contour 
r. Smearing the branches r+ and r_ on the real axis 
and taking into account Eq. (|C1() . the r.h.s. of Eq. (|C3() 
yields the following contribution to the current 

4e7„e-^-*3(/f^/M^j, (C7) 

where the integration over lu has to be understood from 
—oo to +00. Another contribution comes from the first 
term on the r.h.s. of Eq. I|C4|I . By closing the contour 
of the e integration on the complex upper half plane, 
it is non vanishing only if Im[^] > 0. Therefore, only 
the upper branch r_(_ of F contributes. can then be 
smeared on the real axis and one gets 



-4e7Qe 



-ft CJ 



pfiu;f—^}. (C8) 



A similar procedure can be adopted to evaluate the con- 
tribution coming from the second term on the r.h.s. of 
Eq. IIC4|l . One more time we can close the contour of the 
e integration on the complex upper half plane. The pole 
in e = 17 does not contribute since its residue is zero. 
The other pole is in e = C + Ua' and hence one obtains 



-4e7a 3? ■ 



duj f{uj) 



2n uj 



«7 



\ ^ e 

7a 



i(LJ+U^t)t-'yt 



UJ + Ua' - ?7 



(C9) 

Next, we have to calculate the contribution coming from 
Eq. (|C5p . By the same reasoning leading to Eq. ljC9|l it 
is readily verified that it yields the same result. There- 
fore we have to keep in mind that Eq. ljC9|l should be 



multiplied by 2 at the end. Let us now consider the con- 
tribution coming from the first two terms on the r.h.s. of 
Eq. I)C6|I . Since the discontinuous function G'o,o(C) does 
not appear in the integrand we can perform the contour 
integral over C. We find 



4e7o3^ / ^^/(c.'-C/a)- ^. 



- (CIO) 



4e7.K|E7.'/^/0 

I a' 



UJ - Ua') 



UJ + ij 



The contribution coming from the last two terms on the 
r.h.s. ofEq. (|C6I) vanishes. Indeed the integral over e can 
be closed on the complex lower half plane. The pole in 
e — —i"f does not contribute since its residue is zero. The 
other pole contributes only if 3[C] < 0. At the same time 
we can also perform the integration over e' by closing the 
contour in the complex upper half plane. The first term 
in the square brackets of Eq. (|C6|) is non vanishing only 
if 5[C] > 0. The same holds for the second term since the 
pole e' — has vanishing residue. By collecting all the 
results obtained one sees that they can be grouped into 
three broad categories: those which are time independent 
and that give rise to the stationary current, those which 
are proportional to e~^* and those which are proportional 
to e"^'''*. These last ones can be rewritten as 



-4e7.e-^'/^/HE' 



[uj^ + -f^][iuj + Ua,)^ + j^y 

(Cll) 

Let us now group the terms proportional to e Two 
of them comes from Eq. (|C8|) and the first term of Eq. 
(|CT0| : then sum can be written as 



-It 



2tt i(uj + ij)[uj 



(C12) 

The other two pieces come from Eq. (|C9|) (which we 
recall must be multiplied by 2) and the last term of Eq. 
(jClOII . By writing explicitly the real part, after some 
algebra one finds 

^/HE7a'C/a'X (C13) 

Qi' 

wcos[(ct; + C/ct')i] + 7sin[(ti; -I- C/q')^] 

The sum of Eqs. iCTT|l - (|CT2)l - (|CT3| gives exactly the 
quantity Ja{t) ~ jf^ of Eq. 
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